Boundary conditions at the mobility edge 
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It is shown that the universal behavior of the spacing distribution of nearest energy levels at the 
metal-insulator Anderson transition is indeed dependent on the boundary conditions. The spectral 
rigidity E^(_B) also depends on the boundary conditions but this dependence vanishes at high energy 
E. This implies that the multifractal exponent D2 of the participation ratio of wave functions in 
the bulk is not affected by the boundary conditions. 
PACS numbers: 72.15.Rn,73.23.-b,05.45.+b 



The spectral analysis of disordered conductors has 
been proven recently to be a useful tool to probe the na- 
ture of the eigenstates ||l|,||,||,|| . In the diffusive (metallic) 
regime, the conductance g scales linearly with the size L 
of the system, and the wave functions are delocalized over 
the sample. The spectral correlations have been shown to 
be those of random gaussian matrices [g[ , with large de- 
viations above the Thouless energy Ec = h/ro = hD/L^ 
Q , where D is the diffusion coefficient and td is the time 
needed for a wave packet to cross the sample. In partic- 
ular, the distribution P{s) of spacings between nearest 
levels is very well fitted by the Wigner-surmize character- 
istic of chaotic systems 0: P{s) = (7r/2)s exp[— (7r/4)s^] 
where s is written in units of mean level spacing A [pj . 
These deviations of order 1/g'^ M become negligible in 
the limit of large L. In the localized phase, in the limit 
L — > 00, the levels become completely uncorrelated and 
P{s) has the poissonian form: P{s) = exp(— s). This is 
because two levels close in energy are distant in space so 
that their wave functions do not overlap. 

It has been found that the Anderson metal-insulator 
transition in three dimensions is characterized by a third 
distribution ||l| which has the remarkable property to be 
universal, i.e. it is independent of the size, whereas it 
is size dependent in the localized and metallic regimes. 
The transition is thus described as an unstable fixed 
point, in the sense that slightly above the transition 
{W > Wc, where W is the disorder strength and Wc is the 
critical disorder) the distribution tends to a poissonian 
limit when L — > 00, while slightly below the transition 
{W < Wc), it tends to the Wigner-Dyson (WD) distri- 
bution. This third universal distribution has been exten- 
sively studied by several groups who confirmed these re- 
sults, for L ranging from 6 to 100 0,|l|,|l|,|l|,[5|jl6|jl| . 
Up to now, the form of the distribution is still unex- 
plained. 

P{s) carries information on the short range part of 
the spectral correlations. Other characterizations are the 
two-level correlation function (TLCF) of the density of 



states /9(e): K{s) = (p(e + s) p{e)) / {p{e))'^ - 1 and the 
so-called number variance T,'^iE) = {N'^{E)) - (NiE))"^ 
which measures the fluctuation of the number of levels 
N{E) in a band of width E. E is in units of A. It 
is related to the TLCF: E^iE) = 2j^{E - s)K{s)ds. 
Surprisingly enough, the numerical studies which lead to 
the same shape of the distribution P{s) at the transition 
have apparently been all performed using periodic bound- 
ary conditions. In this work, we have calculated P{s) at 
the transition for the same hamiltonian, with different 
boundary conditions (BCs). The hamiltonian is taken as: 



H-E^. 






Cj C^) . 
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The sites i belong to a 3D cubic lattice. Only transfer t 
between nearest neighbours {i,j) is considered. The site 
energies e^ are chosen independently from a symmetric 
box distribution of width W . The metal-insulator transi- 
tion occurs at the center of the band for the critical value 
Wc = 16.5 ± .2 ||I|,|ll|,|l|,0. We have found that although 
the level statistics at the transition is independent of the 
size of the system, it depends on the boundary conditions. 
Our main result is shown in Fig. 1 where we have plotted 
the spacing distribution for four types of BCs, a) periodic 
in the three directions (the situation studied by previous 
authors and that we will refer to as (111), b) periodic in 
two directions and "Hard Wall" (HW) (Dirichlct) in the 
third (110), c) periodic in one direction and HW in the 
two others (100), d) HW in the three directions (000). 
All these distributions are "universal" in the sense that 
they are size independent. The critical point depends at 
most very weakly on the choice of the BCs. It seems 
to shift slightly to smaller W when the number of HW 
boundaries is increased. Using standard scaling analy- 
sis of (s^) and J^ P{s) ds, we found Wc = 16.0 ± .5 for 
the (000) geometry. However, within the range of sizes 
studied (L = 12, . . . , 22) the difference between the P{s) 
at W = 16.0 and at W = 16.5 is negligible compared 
to the remaining statistical fluctuations of the spacing 



distribution. 

In Fig. 2 we have plotted the second moment of the 
level spacing (s^) as a function of the size for the differ- 
ent BCs. This plot shows that the distributions are size 
independent and that they do not converge to a single 
one in the large size limit. 

It may appear a priori surprising that the distribu- 
tion is at the same time size independent and sensitive 
to the BCs. To clarify this point, it is instructive to re- 
call the behavior of the typical dimensionless curvature 
gd = ■"'(IcD/A of the energy levels when an infinitesimal 
flux is introduced in the cylinder geometry. In the metal- 
lic regime, gd{L) increases linearly with the size and it 
decreases exponentially in the localized regime. At the 
transition, the curvature gd{L) — g^ is size independent 
|l8| , [l9| . Since gd measures the sensitivity of energy levels 
to a change of the BCs, the simple fact that it is non-zero 
shows that the spectral correlations can be at the same 
time size independent and sensitive to the BCs. This uni- 
versal sensitivity to the BC has already been discussed in 
the case of periodic BCs where one or several Aharonov- 
Bohm ( AB) fluxes were applied [^|Tj . However in that 
case, the symmetry — time reversal invariance — was at 
the same time broken by the fluxes, such that it is not 
surprising that the statistics is changed. 

The distribution found by other authors with periodic 
BCs is the most rigid of the four distributions that we 
have studied. When periodic BCs are relaxed and re- 
placed by hard wall BCs, the distribution becomes closer 
to the Poisson distribution, with a short range repulsion 
which is characterized by a larger slope of P{s). The 
slope -P'(O) varies by more than a factor three from 2.14 
§ for the BC (111) to 6.80 for the BC (000) (see table 

!)• 

It is useful to stress that in the metallic regime it- 
self, there are deviations to the WD distribution which 
depend on the BCs. These deviations are related to a 
contribution of the diffusive modes 0. At small s, the 
slope of P{s) is given by 



P,s =- 1 + ^^ s, 
6 Tr"^^ 



(2) 



where the coefficient a describes the diffusive motion and 
is given by 






(3) 



For an isolated system, the diffusion modes are quan- 
tized by the BCs. In a direction where the boundaries 
are hard walls, q — nir/L with n = 0, 1,2, 3,---. In a 
direction where the boundaries are periodic, q — 2mr/L 
with n = 0, ±1, ±2, ±3, • • •. In d = 3, one finds: am = 
1.03 , aiio = 2.15 , aioo — 3.39 , aooo = 5.13. So in a 
metal, the slope of P{s) depends on the BC. However, 



the corrections are order l/g"^ and decrease with the size 
since g{L) ~ L, and they vanish for the infinite system. 
At the transition g = g* is size independent and one may 
expect that the correction to P{s) still depends on the 
BC through the quantization of the anomalous diffusion 
modes. This correction can be also simply calculated for 
an anisotropic system. It depends on the shape of the 
sample. This certainly means that the spectral correla- 
tions at the transition are also shape dependent [ p2[ . 

The distributions we have found bear an interesting 
similarity with another recently studied distribution [ p3| . 
A remarquable and simple spectral sequence which is in- 
termediate between the WD and the Poisson distribu- 
tions is obtained by taking the middles of a poissonian 
sequence. This new sequence has been baptized "semi- 
Poisson" [g3|. The corresponding P{s) is given by P3| , p4| : 

-2s 



Pis) = 4se" 



(4) 



It has been shown that the equilibrium distribution 
of charges in a Coulomb gas with logarithmic interac- 
tion only between nearest neighbors is also described by 
eq.(|). The TLCF and T.'^iE) for this model (referred to 
later as Short Range Plasma Model, SRPM) are however 
different from those for the semi-Poisson sequence. We 
shall return to this point later. 

In the inset of Fig. 3 we have plotted the arithmetic 
average of the four distributions. Quite amazingly, it is 
very close to the semi-Poisson distribution. The average 
of the slope at small separation calculated with the four 
BCs is 4.08 ± 0.4 instead of 4 for the semi-Poisson. As 
another characteristic of P{s), the second moment (s^) 
is shown in table I for the various BCs. The average over 
the different BCs is found to be 1.51 ± 0.01. It is 3/2 
for the semi-Poisson. The tails of P{s) have also been 
considerably studied |ll|,|,|jl3|jl|,|l|l . The inset of Fig.l 
shows the tails for the four BCs. They clearly differ by 
the rapidity of their decay, the usual periodic BCs giving 
rise to the fastest decay. It is interesting to notice that 
the behavior at large s is much more affected by HW 
BCs than by a simple addition of a AB flux or even a 
magnetic field po|] . 

We have also investigated random boundary condi- 
tions, with random hopping terms tij = tji connecting 
opposite sides of the sample. Drawing the tij for each 
disorder realization independently from a box distribu- 
tion centered around zero and with width r, we found a 
continuous family of "universal" critical ensembles, which 
are for finite r distinct from the ones with "determinis- 
tic" boundary conditions. 

We now turn to the number variance. A linear behav- 
ior at large E, Y?{E)/E -^ x defines the level compress- 
ibility x^ which is also related to the f ^ dependence 
of the form factor K{t), the Fourier transform of K{s). 
One has X — J^ K{s)ds = K{0) ■ This means x = 1 for 
the Poisson and semi-Poisson sequences, x = for the 
WD correlations and x = 1/2 for the SRPM. 



In Fig. 4, we have plotted Tp'{E)/E for the various 
BCs. It is seen that, hke for P(s), the rigidity depends on 
the BCs for small energy ranges. The rigidity is weaker 
for non-periodic BCs. However, when E increases, the 
different rigidities seem to converge towards the same 
value (see inset of fig.4). We find x — 0.27 ± 0.02 in 
agreement with previous authors ||iq,M. Within error 
bars, this asymptotic value does not depend on the BCs. 

P{s) and X carry information on different time scales 
in the problem. Let us remind that the metallic spec- 
trum is characterized by two characteristic time scales, 
the Thouless time td and the Heisenberg time th — h/A, 
with th/td = 2TrEc/A ^ 1. At the transition, these 
two times are of the same order. Consequently, corre- 
lation functions like P{s) which probe correlations at 
energy scales of the order of A, i.e. time scales of the 
order of th — tjj, probe the sensitivity to the bound- 
ary conditions of a wave packet evolution. However, the 
asymptotic form of the spectral rigidity at energy E ^ A 
typically probes time scales t <^ tn — tjj for which the 
diffusion of a wave packet is insensitive to the BCs. 

More precisely, the form factor has been shown to be 
related to the return probability P{t) for a wave packet 
p6| , p7| . At the transition, the wave functions have a mul- 
tifractal structure pq,E9[ , with a long range tail showing 
a power law decay. Multifractality is characterized by 
the time dependence, when t <^td'- 



p(t) oc i-■°^/^ 



(5) 



where D2 is the multifractal exponent of the inverse par- 
ticipation ratio, X]r(IV'n(r)|'^) oc L^^'^ |^. From the 
limit t ^ of P{t), X is found to be p^: 



1 



X 



(1 



D^, 



)■ 



The multifractal exponent D2 defined from small time 
behavior is thus expected to be independent of the BCs. 
Therefore x should not depend on the BCs either. This 
is seen to be true from the inset of fig. 4, where we show 
that Y?{E)/E converges to the same value x — 0.27±0.02 
for periodic and HW BCs. This value of x leads to a 
multifractal exponent D2 — 1.4 ± 0.2 which has to be 
compared with the value D2 — 1.6 ± .3 found from other 
direct numerical calculations involving the study of the 
wave functions pO]. 

In conclusion, wc have shown that the spectral cor- 
relations at the metal-insulator transition, although be- 
ing universal, i.e. independent of the size, strongly de- 
pend on the choice of the boundary conditions. This 
dependence is most pronounced for small energy scales. 
When the periodic BCs are replaced by hard walls in one 
or more directions, the spectrum becomes less and less 
rigid. However, the level compressibility defined from the 
E —> 00 limit of the spectral rigidity is independent of 
the choice of the boundaries. 
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Table I : Numerical results for the various measures 
of spectral correlations compared with the SRPM. The 
relative errors (standard deviations from 6 system sizes, 
L = 12, 14, 16, 18, 20, 22; 500 to 33 disorder realisations) 
for B are always less than 10% and for s"^ less than 1%. 
(^) This is the value deduced from the Wigner surmize; 
(2) See also ref. §; (3) See also ref. [|l|. 




FIG. 1. Distribution P{s) at the metal- insulator transition 
with four diflterent types of boundary conditions defined in 
the text: A 111, o 110, □ 100, and o 000. Distributions 
with L = 8 to L = 14 are shown. The Wigner-Dyson result 
(continous line) is also plotted. In the inset the tails of P{s) 
are shown for L = 10 and compared with the Semi-Poisson 
distribution, eq.(4) (dashed line). 
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FIG. 4. Y?{E)IE for the diflterent BCs (average over 
1/ = 20 and 22). Symbols as in Fig.l. The inset shows that at 
large energy the difference between the BCs (111) and (000) 
vanishes. 
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FIG. 2. (s ) versus linear size L for different BCs (symbols 
as in Fig.l) for W = 16.5. 




FIG. 3. P{_s) for the four different BCs compared with 
semi-Poisson (dashed line). Inset: "Average" P(s) at the tran- 
sition (v)i compared with semi-Poisson. 
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